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Based on ah initio density functional theory, we have investigated a microscopic mechanism that 
leads to Fcg^/Fcg^ dimerization and consequent ferroelectricity in charge ordered Fe3 04 with P2 
symmetry. In addition to the simple inter-site Coulomb repulsion, quantum hybridization effects 
are invoked to explain the Feg^/Fcg^ bond dimerization. Our results, based on the generalized 
gradient approximation + Hubbard U (GGA+?7) method, indicate that noncentrosymmetric P2 
magnetite shows a finite and sizeable ferroelectric polarization along the h crystalline axis. From 
the U dependence of polarization, we conclude that the origin of ferroelectricity in P2 Fe3 04 lies in 
the recently proposed "intermediate site/bond-centered charge ordering". 

PACS numbers: 75.50.Gg, 77.80.-e 



I. INTRODUCTION 

Advanced multifunctional materials showing coexist- 
ing and spontaneous long-range dipolar and magnetic 
orders are nowadays referred to as multiferroics [1]; they 
combine extremely fascinating physics phenomena with 
a huge potential for technological applications. Recently, 
large interests have been devoted to improper multifer- 
roics [H, [H: here, ferroelectricity is driven by unconven- 
tional mechanisms, such as the loss of centrosymmetry 
due to the occurrence of long-range spin or charge or- 
der (at variance with more conventional ferroelectrics in 
which the ionic degrees of freedom play a prominent role 
in inducing permanent dipole moments). 

Magnetite, Fe304 (formally Fe^+pe^-^+Fe^-^+jBO^J), 
crystallizes in the inverted cubic spinel structure Fd3m 
and shows a metallic behavior at room temperature. Fca 
and Fcb ion sites are coordinated to O ions in the inverted 
spinel structure Fe304, i.e., the tetrahedrally surrounded 
Fe sites are occupied by Fca atoms, whereas the octahe- 
drally coordinated Fe sites are occupied by Fee atoms. 
The latter form a pyrochlore lattice, which consists of 
a network of corner sharing tetrahedra. Ferrimagnetism 
is stable in Fe304, with magnetic moments of Fca sites 
antiparallel to those of Fee sites. Fe304 undergoes a 
first order metal-insulator transition (called Verwey tran- 
sition) at around 120 K, where the resistivity of Fe304 in- 
creases by two orders of magnitude. Q This phase tran- 
sition accompanies a change in the crystal structure of 
Fe304, from cubic to monoclinic symmetry. Verwey pro- 
posed that the metal-insulator transition originates from 
a charge ordering at the Fee sites (Fe^^ [Fe^+Fe^+]B04~), 
similar to strongly correlated hole-doped manganites. 
[1] A full charge localization with ideal Fe|^ and Fe^^ 
does not occur, since the charge disproportionation is 



* Electronic address: |silvia.picQzzi@aquila.infn1t] 



much smaller than the ideal 3+/4+ ionic picture, due 
to strong covalency effects. Anderson pointed out that, 
when putting two Fe^+ and two Fe^+ on each Fee tetra- 
hedron, the number of Fe^+ and Fe^+ pairs is maximized 
and this gives rise to the lowest energy from the Coulomb 
repulsion point of view. However, the Anderson cri- 
terion for low temperature Fe3 04 are inconsistent with 
recent experiment results. [H, 0] As such, the charge or- 
dering occurring at the Verwey transition in Fe304 is still 
under discussion. 

Ferroelectricity has been recently observed at low tem- 
peratures in Fe304 and attracted great scientific inter- 
ests, [l, [il Khomskii et al. have proposed that Fe304 
might show ferroelectricity by means of an intermediate 
site/bond-centered charge ordering, being one of the first 
multiferroic materials. [10] We investigated ferroelectric- 
ity in Fe304 by ah initio calculations, as showing that 
magnetite has non-centrosymmetric Cc symmetry struc- 
ture in the ground state and the origin of ferroelectricity 
in Fe304 is not the site/bond-centered charge ordering, 
but a "charge shift" between specific Fcb site. [ll| 

In contrast to our previous study, in this paper we 
focus on noncentrosymmetric P2 (No. 3) structure in or- 
der to realize the picture proposed by Khomskii, where 
the ferroelectric polarization arises from the Fe|+/Fe|+ 
dimerization. We elucidate three physical properties in 
Fe304 based on ah initio density functional theory. The 
first point is a mechanism for Fcg^/Fcg^ dimerization 
along the monoclinic h direction. Our results indicate 
that the hybridization of electronic orbitals between Fe|^ 
sites {i.e. a quantum mechanical effect) is as important 
as the classical inter-site Coulomb interactions for the 
Fe|^/Feg^ bond dimerization. The second point is re- 
lated to ferroelectricity in P2 Fe304 as induced by in- 
termediate site/bond-centered charge ordering proposed 
by Khomskii et. al. [TO] Ferroelectric polarizations (cal- 
culated via the Berry-phase approach) are compared to 
those obtained in Ref.[ll| for Cc magnetite. Finally, we 
investigate the stability of antiferroelectric P2/c (No. 13) 
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FIG. 1: Illustration of a net polarization (red arrow) by the 
intermediate site-centered and bond-centered charge ordering 
in Fes 04- The Fe^^ and Fe^^ charge ordering and an alter- 
nation of short and long bond lengths between Fe^^ and Fe^^ 
sites in low temperature magnetite are shown. 



and ferroelectric P2 Fe3 04 vs electronic correlation ef- 
fects. 




FIG. 2: Schematics of the corner sharing tetrahedra formed by 
Fee sites in P2/c Fe3 04. The labels of Fes site are consistent 
with Refill. Black and white squares indicates Fe|^ and Fe^^ 
sites, respectively. 



II. ELECTRIC POLARIZATION INDUCED BY 
CHARGE ORDERING: BASIC CONCEPTS 

Khomskii et al. have discussed how charge ordering 
could induce ferroelectricity in magnetic systems. [1^ 
According to their idea, when either site-centered or 
bond-centered charge orderings happen individually, a 
system does not have a net dipole moment: however, if 
the site-centered and bond-centered charge orderings oc- 
cur simultaneously in a system, the inversion symmetry 
is broken so that the system shows a net dipole moment, 
therefore becoming ferroelectric {of. FigH]). [10] Based 
on this scenario, Khomskii et al. have proposed that 
Fe304 could become ferroelectric, with the monoclinic 
b axis as polarization direction. Actually, experimen- 
tal results indicate that the situation in Fe304 is similar 
to the intermediate site/bond-centered charge ordering: 
Wright et al. [12] have investigated the crystal structure 
of Fe304 at low temperature by using high-resolution 
neutron powder-diffraction and X-ray powder-diffraction 
measurements. According to Ref.[l2|, the distances be- 
tween charge ordered Fe|^/Feg^ sites along b are strongly 
modulated, from 2.86 A to 3.05 A. So, in addition 
to the site-centered charge ordering of Fe^^ and Fe^^ 
sites, there is a sequence of short and long FeB-Fee bond 
lengths along b {i.e. bond dimerization process). This 
situation indeed corresponds to the site/bond-centered 
charge ordering proposed by Khomskii et al [l^., the so 
called "Zener polaron" . 

Generally, the origin of bond dimerization in other 
multiferroic materials is the spin configuration, e.g., in 
HoMn205 the bond lengths between antiparallel (paral- 
lel) Mn ions are expanded (contracted) so as to optimize 
the double exchange mechanism. [13] There is, however, 
no correlation between the dimerization process and the 
spin configuration in Fe304, since all the spin directions 
of Fcb sites are the same. Wright et al. have proposed the 
bond length modulation to be due to the simple classi- 



cal inter-site Coulomb repulsion between charge ordered 
Fcg^ and Fe^^ sites, i.e., below the Verwey tempera- 
ture the corner sharing Fee tetrahedra are constructed 
by three Fe^^ and one Fe|^ ions or by three Fe|^ and 
one Fcg^ ions (3-1 charge order), at variance with the An- 
derson criterion in which the configurations of two Fe^^ 
and two Fe^^ ions (2-2 charge order) are most stable [5[, 
as shown in Fig|2l It can be easily understood that a 
system can lower its energy, from the inter-site Coulomb 
interaction point of view, if Fee ions on B(3) and B(4) 
sites move along the thin arrows in FigO As a result, 
Fcb ions on B(3) and B(4) sites are moved along the bold 
arrows in FigO the net effect being a bond dimerization 
between B(3) and B(4) sites. However, as we'll show, an 
important role in the dimerization process in Fe304 is 
played by the on-site Coulomb interaction of Fe atoms 
as well, and understanding the dimerization process in 
Fe304 is not so trivial. 



III. CALCULATION METHOD 

Electronic structure calculations and structural op- 
timizations are performed by the "Vienna Ah initio 
Simulation Package" (VASP) and Projector Augmented 
Wave (PAW) pseudopotentials. |1^] The Perdew-Becke- 
Erzenhof (PBE) of the generalized gradient approxima- 
tion (GGA) is employed for the exchange-correlation po- 
tential. [l5[ The plane wave cutoff energies are 400 eV for 
Fe and O atoms. The GGA+/7 calculations within Du- 
darev's approach are performed by_applying a Hubbard- 
like potential for Fe-d states only, [l^ Effective Coulomb 
energies U = 4.5, 6.0, and 8.0 eV and an exchange pa- 
rameter J = 0.89 eV are used. [13] To calculate the 
electronic structure of ferroelectric Fe304, we consider 
an "artificial" crystal structure with P2 symmetry, built 
by copying the lower half of the (experimentally sug- 
gested) monoclinic P2/c antiferroelectric cell to the up- 
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FIG. 3: Red arrows show local polarization directions in (a) 
antiferroelectric P2/c and (b) ferroelectric P2 Fes 04- Purple 
and orange spheres indicate Fe^^ and Fe^^ ions, respectively. 
For clarity, Fca and O ions are not shown. Tetrahedra com- 
posed by three Fe^^ and one Fe^^ ions and three Fe^^ and 
one Fcg^ ions are highlighted by orange and purple planes, 
respectively. 



per half, as shown in FiglSl Note that there is no net 
polarization in antiferroelectric P2/c magnetite, because 
the local polarizations in the upper half and lower half 
cancel each other. On the other hand, ferroelectric P2 
Fe304 can show a finite net polarization. Computa- 
tional details for the electronic structure and structural 
optimization in ferroelectric P2 Fe304 are the same as 
antiferroelectric P2/c. Both P2/c and P2 contain 56 
atoms (24 Fe and 32 O ions) in a^/ V^xad V^>^2ac unit 
cell, where stands for the lattice constant of cubic 
Fd3m Fe304. Lattice parameters are fixed to experi- 
mental values: a = 5.9341 A, 6 = 5.9256 A, c = 16.7524 
A, and (3 = 90.2000°. [IH In the structural optimiza- 
tion, a threshold on the atomic forces of 0.01 eV/A is 
employed. Internal atomic coordinates are optimized in 
P2/c and P2 Fe304 starting from experimental Wyck- 
off parameters. |12| Fe-3p^3d^4:S^ and 0-2s^2p^ electrons 
are treated as valence electrons. For Fe atoms, the ferri- 



magnetic configuration is considered: all Fca sites have 
up-spin and all Fcb sites have down-spin. We neglected 
spin-orbit coupling, so that the direction of magnetiza- 
tion with respect to the crystal is not specified. For the 
P2/c and P2 Fe304, the 6x6x2 Monkhost-Pack fe-point 
grid in the Brillouin zone is used. The Berry phase ap- 
proach [l^, [iO] within the PAW formalism and the point 
charge model with nominal charges (Fca: 3+, Fe|^: 2-h, 
Fcg^ : 3+, and O: 2-) are employed to calculate polar- 
ization. We integrate over six /c-point strings parallel to 
the polarization direction (crystallographic b direction) 
in the Berry phase approach. 

IV. RESULTS 

A. Fe^+/Fe^+ bond dimerization in Fe304 

First, we will show that the site/bond-centered charge 
ordering occurs in Fe304, so that ferroelectricity can be 
expected. Figure [H illustrates charge densities (in a se- 
lected energy range, see caption) of Fe|^, Fe^^, and 
sites in ferroelectric P2 Fe304 and an enlarged view of 
one of the z=l/8, 3/8, 5/8, and 7/8 planes in FiglSKb). 
As discussed in a previous study ^ , one electron occu- 
pies one of the spin down t2g orbitals {xy^ yz^ and zx) 
on Fe^+ sites {d^ configuration, ^2^T^^T^2c/i)5 therefore, 
charge densities on Fe^+ sites have a cross-like shape 
corresponding to t2g orbitals and extend to interstitial 
regions of anions 0^~. The spin down t2g orbitals on 
Fe^+ sites {d^ configuration, t^g^e'^^) are pushed above 
the Fermi level. The important point here is the charge 
density on the sites: charge is not present on 
sites in the long bonds (c/r. FigH]), whereas in the short 
bond sides there is an appreciable charge density with 
Oxygen p character. The latter can be viewed as "bond- 
centered charge". It can be easily seen that there is no 
overlap between d orbitals on Fe|^ and Fe^^ sites and 
p orbitals on sites in the long bond sides; however, 
the d orbitals between Fe|^ and Fe^^ sites via p orbitals 
on sites can hybridize on the short bond sides. Note 
that the hybridization between d orbitals on Fe|^ sites 
and p orbitals on sites is not perfectly 7r-like because 
of the small tilt of p orbitals on sites, so that we can 
expect a somewhat larger hybridization with respect to 
the ideally undistorted atomic arrangement. As a result, 
due to orbital hybridization, short bond sides have higher 
electron density and long bond sides have lower electron 
density. The center of gravity of charges is therefore ex- 
pected to move away from the ions and we infer from 
FigHlthat, in addition to the 5i^e-centered charge order- 
ing on Fcb sites, even the &on(i-centered charge ordering 
proposed by Khomskii et al. is realized in Fe304. [T^ 

Table [J shows the dependence of Fee -Fee bond lengths 
on the on-site Coulomb interaction U in the corner shar- 
ing tetrahedra and valence charges qd of Fee sites in fer- 
roelectric P2 Fe304 (the numbering of Fcb sites corre- 
sponds to that of P2/c structure in Figl2|). As shown in 
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FIG. 4: Minority-spin charge densities of Fe^ , Fe^ , and 
sites in ferroelectric P2 Fe3 04. A closeup of one of the 
z=l/8, 3/8, 5/8, and 7/8 planes in Fig|3];b). Short and long 
mean alternations of bond lengths between Fes sites due to 
the Fe^^/Fe^^ bond dimerization along b. Charge densities 
are calculated within 1.00 eV below the Fermi level. 
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FIG. 5: Schematics of "asymmetric hybridization" between 
Fcg^ sites in the minority spin state. Thin black arrows de- 
note the electron spin and thick blue arrows indicate electron 
hopping. 



the upper part of Table [J and FiglSl there are alterna- 
tions of short and long bond lengths between Fee sites, 
leading to the Fe^+/Fe^+ bond dimerization along the b 
axis, so that the corner sharing tetrahedra formed by Fee 
sites are strongly distorted from ideal regular tetrahedra. 
Table |I] also indicates that the distortions of the tetra- 
hedra are suppressed upon increasing U. For example, 
along the b direction, the bond lengths between FeB(3) 
and Fee (4) sites are modulated from 2.888 to 3.038 A 
for [7=4.0 eV, from 2.919 to 3.006 A for U=6.0 eV, and 
from 2.934 to 2.992 A for [7=8.0 eV, respectively (the 
bond length is modulated from 2.880 to 3.046 A in the 
experiment [7]). 

Now, we focus on a mechanism of the Fe^^/Fe^^ bond 
dimerization in Fe304. Previous studies [l2|, [13] sug- 
gested dimerization on Fe|^ and Fe^^ sites to origi- 
nate from the inter-site Coulomb interaction between Fee 
sites. However, this is not the only factor and we pro- 
pose that "asymmetric hybridization" plays a relevant 
role in this context (c/r. Fig. [5]). As mentioned above, 
one electron occupies the t2g down spin states at Fe^^ 



sites, whereas t2g and Cg down spin states at Fe^ sites 
are located above the Fermi level. In such conditions, 
Fe|^ electron orbitals can hybridize, so that the Fe|^- 
Fe|^ bond lengths are shrinked. However, electron or- 
bitals on Fe|^ and Fe^^ sites hybridize much less, due to 
the larger energy level difference between the t2g down 
spin states in Fe|^ and Fe^^ sites; therefore, shortening 
of bond lengths between Fe|^ and Fe^^ sites can not be 
expected by orbital hybridization. In this sense, we refer 
to the hybridization of electron orbitals on Fe|^ sites in 
Fe304 as "asymmetric hybridization". 

In order to prove the validity of the "asymmetric hy- 
bridization" , we investigate the on-site Coulomb interac- 
tion dependence of valence charges qd of Fcb ions (calcu- 
lated within a Wigner-Seitz radius r = 1.00 A). As shown 
in the lower part of Table [H the Fe|^/Feg^ charge or- 
dering occurs in Fe304. This Fcp^/Fe ^ c harge ordering 
has been well known and investigated. [13] We note that, 
when U is set to zero, there is neither charge order nor 
bond dimerization of Fee sites, due to metallicity. The 
key point here is the U dependence of charge separations 
between Fcb sites. Our results indicate that, upon in- 
creasing [7, the charge separation in Fee sites increases, 
due to the fact that electrons are more localized for large 
U. Therefore, upon increasing the charge separations 
in Fcb sites, one expects the inter-site Coulomb inter- 
actions to more strongly distort the Fcb tetrahedra, im- 
plying a larger tendency towards bond dimerization be- 
tween Fcb sites. However, Table |T] shows that the dimer- 
ization on Fee sites decreases when increasing U. We 
infer that these results can not be explained by inter-site 
Coulomb interaction only. Indeed, from the "asymmet- 
ric hybridization" point of view, it can be expected that, 
due to the stronger localization of electrons for large [7, 
the reduced hybridization leads to a weakening of dimer- 
ization on Fcb sites (as confirmed by Table H]) . 



B. Ferroelectricity in Fe304 

From the discussion in the previous section, we can 
expect ferroelectric polarization along b in P2 Fe304 to 
occur by means of the site/bond-centered charge pro- 
posed by Khomskii et. al. jl3] Due to the presence of 
inversion symmetry, P2/c Fe304 doesn't have a net po- 
larization; however, there is no inversion symmetry in 
P2 Fe304, so the system is allowed to have finite polar- 
ization. We choose the antiferroelectric P2/c Fe304 as 
a centrosymmetric reference structure, i.e., spontaneous 
polarization is obtained by polarization difference be- 
tween ferroelectric P2 and antiferroelectric P2/c Fe304, 
AP = PpE — ^AFE- Here, polarization along b in Fe304 
is calculated by two different methods: one is the Berry 
phase method, where the macroscopic polarization is ob- 
tained from the Wannier- center displacements, [19] an- 
other is the simple "point charge model" with nominal 
charges on the Fe an O ions. Note that the polar iza- 
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U (eV) B(la)-B(lb) B(la)-B(3) B(la)-B(4) B(lb)-B(3) B(lb)-B(4) B(2a)-B(2b) 



4.5 


2.9671 


2.9215 


2.9129 


2.9118 


2.9046 


2.9671 


6.0 


2.9671 


2.9169 


2.9366 


2.9067 


2.9289 


2.9671 


8.0 


2.9671 


2.9156 


2.9510 


2.9058 


2.9431 


2.9671 


U (eV) 


B(2a)-B(3) 


B(2a)-B(4) B(2b)-B(3) B(2b)-B(4) 


B(3)-B(4) 


B(4)-B(3) 


4.5 


3.0028 


3.0195 


3.0091 


3.0247 


3.0379 


2.8877 


6.0 


2.9958 


3.0080 


3.0017 


3.0129 


3.0064 


2.9193 


8.0 


2.9920 


2.9977 


2.9984 


3.0031 


2.9923 


2.9336 


U (eV) 


B(la) 


B(lb) 


B(2a) 


B(2b) 


B(3) 


B(4) 


4.5 


5.373 


5.372 


5.140 


5.140 


5.166 


5.358 


6.0 


5.398 


5.397 


5.106 


5.106 


5.116 


5.401 


8.0 


5.413 


5.412 


5.073 


5.073 


5.077 


5.416 



TABLE I: Upper and middle parts show the on-site Coulomb interaction U dependences of FeB-Fee bond lengths (in A) in 
the corner sharing tetrahedra. The lower part shows valence charges qd (in e) of Fee sites in P2 Fe304. In the middle table, 
B(3)-B(4) and B(4)-B(3) indicate long and short bond sides, respectively. 
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^7=4.5 eV U=6.0 eV U=S.O eV 



FIG. 6: Dependence of electric polarization in ferroelectric 
P2 Fe304 on the on-site Coulomb interaction. Considered U 
values are 4.5, 6.0, and 8.0 eV. The exchange parameter J 
is fixed to 0.89 eV. Polarizations calculated by point charge 
model and by the Berry phase method are shown. 



centers. Such enhancement of polarization is of course 
not taken into consideration in the point charge model. 
Additionally, calculated polarizations decrease when in- 
creasing the U value. The obtained results are consis- 
tent with the Feg^/Fcg^ bond dimerizations discussed in 
Sec lIV A[ As increasing [/, the Fe|^/Feg^ dimerizations 
are suppressed, so it can be easily understood that if 
the bond dimerization between Fe|^ and Fe^^ sites dis- 
appears, polarization caused by the site/bond-centered 
charge ordering vanishes. Another important point here 
is the polarization with large U values. As shown in 
Figini polarizations calculated by the point charge model 
and the Berry phase method for large U approach the 
same value. This can be understood as follows: upon 
increasing /7, electrons are more localized on Fe ions. 
Therefore, since the Wannier centers get closer to Fee 
sites and there is no enhancement of dipole moments 
by the displacements of the Wannier centers, the point 
charge model and the Berry phase method have similar 
values for large U. 



tion can only be determined up to an integer multiple 
of the polarization quantum |23|]. In ferroelectric P2 
Fe304 structure, the polarization quantum eR/Vt is 16.12 
/iC/cm^ along h direction, where e is the electron charge 
and R is the lattice vector along h. 

Figure [6] shows polarizations along h obtained by point 
charge model and by the Berry phase method. In the 
later approach, the internal atomic positions of ferroelec- 
tric P2 and antiferroelectric P2/c Fe304 are optimized 
in the conditions of U =4.5, 6.0, and 8.0 eV, and then the 
same U values as for the crystal optimizations are consis- 
tently used for polarization calculations. At /7=4.5 eV, 
the polarization calculated by the Berry phase method is 
much larger than that of the point charge model. This 
large difference is due to strong enhancement of dipole 
moments by the Wannier centers displaced from atomic 



We now comment on the difference between our present 
results and our previous calculations. [ll| We showed 
that Fe304 has Cc symmetry in the ground state and 
polarization direction is not along h but in the monoclinic 
ac mirror plane, PBerry=(-4.41, 0, 4.12) in /iC/cm^ (for 
[7=4.5 eV and J=0.89 eV). The difference between our 
present results and previous calculations originates from 
the microscopic mechanism of ferroelectricity: in Fe304 
with Cc symmetry, the polarization is not driven by the 
site/bond-centered charge ordering but by "charge shift" 
at selected Fcb sites. The bond dimerization on Fe^^ and 
Fcg^ sites along 6 in Cc Fe304 does indeed exist; however, 
dipole moments caused by the site/bond-centered charge 
ordering are cancelled between upper half and lower half 
of the unit cell, due to the Cc symmetry. Therefore, the 
site/bond-centered charge ordering is not appropriate as 
a mechanism of ferroelectricity in Cc magnetite. 
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FIG. 7: On-site Coulomb interaction dependence of total en- 
ergy difference per unit cell between antiferroelectric P2/c 
(AFE) and ferroelectric P2 (FE) Fe304. 



C. Stability of antiferroelectric P2/c and 
ferroelectric P2 Fe304 

Finally, we briefly discuss the stability of antiferroelec- 
tric P2/c and ferroelectric P2 Fe3 04. Figure [3 shows the 
total energy differences per unit cell of antiferroelectric 
and ferroelectric Fe304 (after structural optimizations) 
as a function of on-site Coulomb interaction U. The 
structural optimizations are performed separately for fer- 
roelectric and antiferroelectric Fe304 for [7=4.5, 6.0. 8.0, 
and 10.0 eV. Incidentally, we note here that at [7=4.5 
eV, Fe304 with Cc symmetry is the ground state pj| . 
When the on-site Coulomb interaction is set to U = 4.5 
eV, the ferroelectric state is less stable than antiferro- 
electric state as shown in Fig|71 However, the situation 
is changed for large U values. The energy difference be- 
tween antiferroelectric P2/c and ferroelectric P2 Fe304 



decreases with increasing U value; when U reaches ~8.0 
eV the ferroelectric state becomes even more stable than 
the antiferroelectric state, paving the way to P2 mag- 
netite as a possible multiferroic. 



V. SUMMARY 

Microscopic mechanisms leading to Fe|^/Feg^ dimer- 
ization and ferroelectricity in Fe304 have been investi- 
gated by ab initio density functional theory with Hub- 
bard U (GGA+/7) method. From the on-site Coulomb 
interaction U dependence of Fcb-Fcb bond lengths in 
corner sharing tetrahedra, we show that the "asymmet- 
ric hybridization" between Fe|^ sites plays an impor- 
tant role in Fe^ /Fe^ bond dimerization, in addition to 
the simple inter-site Coulomb interaction. Additionally, 
based on the results of the U dependence of polariza- 
tions and Fcb-Fcb bond lengths, we show that the non- 
centrosymmetric Fe304 with P2 symmetry can become 
ferroelectric by means of the intermediate site/bond cen- 
tered charge ordering proposed by Khomskii, where the 
Fe|^/Feg^ site charge order and Fe^ /Fe^ bond dimer- 
ization along the monoclinic b direction simultaneously 
occurs. P2 Fe304 can therefore be well considered as 
an improper multiferroic material with charge-order- 
induced ferroelectricity. 



VI. ACKNOWLEDGMENTS 

The research leading to these results has received 
funding from the European Research Council under 
the European Community's 7th Framework Programme 
(FP7/2007-2013) / ERG grant agreement n. 203523. We 
thank Prof. Daniel Khomskii for useful discussions. The 
crystal structures and the charge states are plotted by 
using the software VESTA 0|. 



[1] For comprehensive reviews, see S.W.Cheong and 
M.Mostovoy, Nature Mater.6, 13 (2007); R. Ramesh and 
N.A.Spaldin, Nature Mater. 6, 21 (2007) and references 
therein. 

[2] S. Picozzi and C. Ederer, J. Phys.: Condens Matter 21, 
303201 (2009). 

[3] E. J. W. Verwey Nature (London) 144, 327 (1939). 

[4] E. J. W. Verwey and P. W. Haayman, Physica (Amster- 
dam) 9, 979 (1941). 

[5] P W. Anderson, Phys. Rev. 102, 1008 (1956). 

[6] J. Garcia and G. Subias, J. Phys. Condens. Matter 16, 
R145 (2004). 

[7] J. P. Wright et. al, Phys. Rev. Lett. 87, 266401 (2001). 
[8] C. Medrano, M. Schlenker, J. Baruchel, J. Espreso, and 

Y. Miyammoto, Phys. Rev. B 59, 1185 (1999). 
[9] M. Alexe, M. Ziese, D. Hesse, P. Esquinazi, K. Yamauchi, 

T. Fukushima, S. Picozzi, and U. Gosele, Adv. Mater. 21, 



4452 (2009). 

[10] D. Khomskii, Bull. Am. Phys. Soc. A13.00006 (2007); J. 

van den Brink and D. 1. Khomskii, J. Phys.: Condens. 

Matter 20, 434217 (2008). 
[11] K. Yamauchi, T. Fukushima, and S. Picozzi, Phys. Rev. 

B 79, 212404 (2009). 
[12] J. P. Wright et al, Phys. Rev. B 66, 214422 (2002). 
[13] G. Giovannetti and J. van den Brink, Phys. Rev. Lett. 

100, 227603 (2008). 
[14] G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 

(1996). 

[15] J. P. Perdew et al, Phys. Rev. Lett. 77, 3865 (1996). 
[16] S. L. Dudarev et al, Phys. Rev. B 57, 1505 (1998). 
[17] H.-T. Jeng et al, Phys. Rev. B 74, 195115 (2006). 
[18] Zuo et al, Phys. Rev. B 42, 8541 (1990). 
[19] R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 
1651 (1993). 



7 



[20] R. Resta, Rev. Mod. Phys. 66, 899 (1994). 
[21] H.-T. Jeng et al, Phys. Rev. Lett. 93, 156403 (2004). 
[22] I. Leonov, A. N. Yaresko, V. N. Antonov, M. A. Korot in, 
and V. I. Anisimov, Phys. Rev. Lett. 93, 164404 (2004). 



[23] J. B. Neaton et al, Phys. Rev. B 71, 014113 (2005). 
[24] K. Momma and F. Izumi, J. Appl. Crystallogr. 41, 653 
(2008). 



